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Abstract 

Twisting transformations for the heat operator are introduced. They 
are used, at the same time, to superimpose a la Darboux N solitons to a 
generic smooth, decaying at infinity, potential and to generate the corre- 
sponding Jost solutions. These twisting operators are also used to study 
the existence of the related extended resolvent. Existence and uniqueness 
of the extended resolvent in the case of N solitons with TV "ingoing" rays 
and one "outgoing" ray is studied in details. 



1 Introduction 

The Kadomtsev-Petviashvili equation in its version called KPII 

(u t - 6uu Xl + u XlXlXl ) Xl = -3u X2X2 , (1.1) 

is a (2+l)-dimensional generalization of the celebrated Korteweg-de Vries (KdV) 
equation. As a consequence, the KPII equation admits solutions that behave at 
space infinity like the solutions of the KdV equation. For instance, if u\{t,x\) 
obeys KdV, then u(t,x\,X2) — U\(t,x\ + ^ix% — 3/i 2 i) solves KPII for an arbi- 
trary constant fj, £ M. Thus, it is important to consider solutions of (jl.ip that 
decay at space infinity in all directions with exception of a finite number of 
1-dimcnsional rays with behavior of the type of u\ . Even though KPII has been 
known to be integrable for more than three decades [1,2], its general theory 
(involving such non-decaying solutions) is far from being complete. Thus, the 
Cauchy problem for KPII with rapidly decaying initial data was solved in [3,4] 
by using the Inverse Scattering Transform (1ST) method, based on the spectral 
analysis of the heat operator 

C{x,d x ) = -d X2 +d 2 Xi -u{x), x = (x 1 ,x 2 ), (1.2) 
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that gives the associated linear problem for the KPII equation. The standard 
approach to the spectral theory of the operator (|1.2p is based on integral equa- 
tions for the Jost solution Q(x, k), where k G C denotes the spectral parameter, 
or for ^/(a;,k), Jost solution of the dual operator C d . However, it is known 
that in the case of potentials with one-dimensional asymptotic behavior these 
integral equations are ill-defined. 

In order to overcome these difficulties, a resolvent approach was developed 
in [5-14]. In this approach a space of operators A(q) with kernels A(x,x';q) 
belonging to the space of tempered distributions of variables x,x',q 6 K 2 was 
introduced. Ordinary differential operators C{x,d x ) are imbedded in this space 
as operators with kernels 

L(x, x 1 ; q) = C{x, d x + q)S(x — x'). (1-3) 

Due to the additional dependence on the parameter q these operators are called 
extended operators. In this space a generalization of the resolvent of a differen- 
tial operator, called extended resolvent, was introduced, enabling consideration 
of the spectral theory of operators with nontrivial asymptotic behavior at space 
infinity. 

In [11-14] we considered the nonstationary Schrodingcr and heat operators 
with potentials with only one direction of nondecaying behavior. The first step 
in solving the problem was the embedding of a pure one-dimensional spectral 
theory in the two-dimensional one, building the two-dimensional extended re- 
solvent for an operator with potential u{x) = ui(x 1 ). The second step was 
the dressing of this resolvent by an arbitrary bidimensional perturbation of the 
potential u\. Finally, all mathematical entities appearing in the 1ST theory 
were obtained from this dressed resolvent by a reduction procedure, allowing 
the formulation of the direct and inverse problems. Let us also mention that 
the standard spectral theory for the heat operator in the case of potentials non- 
decaying in one space direction was developed in [15], but under some special 
and inexplicit conditions on the potential. 

Here we consider the substantially more complicated problem: the case of 
a potential u non decaying along multiple non parallel rays. Thus now there 
is no analogy with the one-dimensional case and the whole theory has to be 
constructed directly, without embedding one-dimensional entities in two di- 
mensions. Therefore, we are obliged to consider directly true bidimensional 
potentials as it was already done in [16, 17] for the nonstationary Schrodinger 
operator. In fact, we apply the same procedure that was successfully applied in 
that case. Precisely, in order to get the potential corresponding to N solitons 
"superimposed" to a generic smooth decaying potential and the related Jost so- 
lutions we construct these entities directly by means of twisting transformations 
instead of using recursively the (binary) Darboux transformation. In this way 
we recover potentials and Jost solutions that were obtained in [11] by using a 
recursive procedure and, also, some alternative representations of the pure soli- 
tonic potentials obtained in [18-22] by using the r functions. It is worth noting 
that the heat operator is not self-dual and thus singular behaviors of the left 
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and right twisting operators are not obliged to be correlated. This explains why 
the structure of the N soliton solution of KPII is much richer than for KPI. In 
particular, the N soliton solution can have a different number of "ingoing" and 
"outgoing" (in the sense of xi — > —oo, or +00) rays. 

The article is organized as follows. In Sec. [5] we briefly review some aspects 
and basic ideas of the extended resolvent approach, referring to the example 
of an operator (|1.2[) with a smooth, rapidly decaying potential u(x). For fur- 
ther details, we refer the interested readers to [5-9]. In Sec.|3]we introduce the 
twisting operators, i.e., operators that "twist" the operator L, extension in the 
sense of (|1.3I) of the operator C in (|1.2p , to an operator L' of the same kind 
with a potential u' describing N solitons "superimposed" a la Darboux to the 
background potential u. In Sec. 0] by using these twisting operators, under the 
assumption of the existence of the resolvent M'(q) we derive an explicit expres- 
sion for its kernel M'(x, x';q). So the problem of the existence of the resolvent 
M'{q) is reduced to the problem of finding the region in the g-plane, where this 
kernel is a tempered distribution. This problem resulted to be more arduous 
than in the case of the nonstationary Schrodinger operator and, therefore, in 
Sec.[5j we first consider the pure iV-soliton potentials and, afterwards, in Sec. El 
the special simple subclass of iV-soliton potentials with N "ingoing" rays and 
only one "outgoing" ray. We prove that for N > 1 the resolvent M'(q) exists 
only in the region of the g-plane exterior to a polygon with N + 1 sides. More 
specifically, we prove that there exists a value ko of the spectral parameter k, 
for which the dual Jost solution ^'(x, k) is such that #'(x, k ) exp(qia;i + (72^2) 
is exponentially decaying on the x-plane for any value of q belonging to this 
polygon and, consequently, we deduce that inside the polygon the twisted op- 
erator L'(q) has a left annihilator and cannot have right inverse. Nonetheless, 
the Green's function of C! exists and can be uniquely derived via a reduction 
from M'(q). In a forthcoming publication we plan to consider the generic case 
of an arbitrary number of incoming and outgoing rays and to elucidate the role 
of these annihilators in the spectral theory of such potentials. 



2 Background theory 

Let us introduce the space of extended operators A(q), i.e., operators with 
kernel A(x, x';q) belonging to the space S' of tempered distributions of the six 
real variables x = (xx,X2), x' = (x^x^), and q = ((71,92)- For two extended 
operators A(q) and B(q) with kernels A(x,x';q) and B(x,x';q) we introduce 
the composition law 

(AB)(x, x';q) = J dx" A(x, x"; q) B{x" , x'; q), (2.1) 

provided the integral exists in terms of distributions. An operator A can have an 
inverse, A^ 1 , in the sense of this composition: AA~ X = I and A^ 1 A = I, where 
I is the unity operator, i.e., the operator with kernel I{x, x'; q) — 8{xi—x' 1 )5(x2 — 
x' 2 ). A specific subclass of extended operators is given by the extensions L(q) 
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of differential operators as defined in (| 1 . 3[) . where C{x, d x ) denotes a differential 
operator whose coefficients are smooth functions of x. Let us associate to any 
operator A(q) with kernel A(x,x';q) its "hat"-kernel 

A(x, x'; q) = e qix - x '^A(x, x';q), (2.2) 

where qx = q%x% + (72^2- For a differential operator L(q) this procedure is the 
inverse of the extension introduced in (|1.3p . i.e., L(x,x';q) = L(x,x'). Notice, 
however, that, since the kernels of the operators A are only subjected to the 
requirement of belonging to the space of tempered distributions, in general, the 
hat-kernel A still depends on q and is not necessarily bounded. The conjugate 
A*(q) of an operator A{q) is defined by 



A*(x,x';q)=A(x,x';q). (2.3) 

It is convenient to introduce the representation of the operator A{q) in the 
p-space which is defined by the Fourier transform 



A(p;q) = J dx J dx' e^^*-^*' A{x,x';^), 



(2.4) 



where p — (pi,p2) G R 2 and we introduced a two dimensional complex vector 

q = qsR+*qo> 9 = q^, q^a^R 2 . (2.5) 

The composition (|2.f |) in the p-space becomes a sort of shifted convolution 

(AB)(p-q) = J dp'A(p-p';q+p')B(p';q). (2.6) 



Then, with these notations, the extension of the heat operator (jTHj) reads 

L = Lq — u, u{x, x'; q) — u(x)5(x — x'), (2-7) 

where Lq is the extension in the sense of (|I.3[) of the differential part Cq of the 
heat operator (|f ,2p and has, in the p-space, thanks to (|2.4[) . kernel given by 

Lo(p;q) = (iqa-q?)*(p). (2.8) 

The main object of our approach is the extended resolvent (or resolvent, for 
short) M(q) of the operator L(q), which is the inverse of the operator L in the 
sense of composition (12.11) (or (|2.6p ). i.e., 

LM = ML = I, (2.9) 

and can be defined as the solution of the integral equations 

M = M + Al uM, M = M + uMM , (2.10) 

where Mq is the resolvent of the zero potential (bare) operator Lq. Thanks 
to (|2.8p we have Mo(p;q) = 8(p)/(i<\ 2 —q\) for the kernel of this operator in 



4 



the p-space. In the case of a rapidly decaying potential u(x), the existence and 
uniqueness of the solution of the equations in (|2.10[) can be proved in analogy 
with [4]. Notice that for a real potential u(x] both L and M are self-conjugate 
operators in the sense of definition (|2.3p . 

Thanks to (|2.6[) and the explicit form of Mo(p;q) given above the integral 
equations (12.10|) written in the p-space show that the kernel M(p; q) is singular 
for q 2 = — iqf and for q 2 +P2 = — *(qi+Pi) 2 - Therefore, it is natural to 
introduce the following truncated and reduced values of the resolvent 

f(f>;q) = (ML )(p;q) , uj{p; q) = (L M)(p; q) 

(2.11) 

One can show that the operator L(q) and its resolvent M(q) admit the following 
bilinear representations in terms of v and u> 

L = vL oj, M = vMquj. (2.12) 

Correspondingly, we call the operators v and to dressing operators since they 
"dress" the bare operators L and M . Notice that f(p;q) and ui(p;q) satisfy 
asymptotics liniq^oo v(p\ q) = S(p) and lim q -kx, o;(p; q) = S(p) and do not 
depend on q 2 , which, if necessary, we make clear by writing v(j>\ q-jj and u)(p; qj). 
The dressing operators v and u> satisfy the equations 

Lv = vLq, u>L = Lquj, (2-13) 

and are mutually inverse 

wv = I, vuj = I. (2.14) 

In order to define the Jost solutions by means of the dressing operators we 
introduce the following Fourier transforms: 

X (x,c ll ) = J ' dpe-^vfath), Z{x,q l )= j dpe-^uip-^-px). (2.15) 

Then the Jost and the dual Jost solutions can be defined, respectively, as 

$(x,k) = e- lkxi - ]i2x2 x{x,k), *(x,k) =e lkxi+k2;r2 e(x,k), (2.16) 

where we denoted q : as k in order to meet the standard notations for the 
spectral parameter . Thanks to (|2.13p , they obey the heat equation and its dual 

(~d X2 +d 2 Xi -u(x)Mx,k) = 0, (d X2 +d 2 Xi -u(x))^(x,k)=0. (2.17) 

3 Darboux transformations via twisting opera- 
tors ( and 77 

In order to build a two-dimensional potential describing N solitons superim- 
posed to a generic smooth background, we bypass the recursive procedure used 
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in [11] and, by using the operator formulation introduced in the previous sec- 
tion, we construct directly the final Darboux transformation. The main tool in 
accomplishing this result is what we call the twisting operators. 

Let us consider a transformation from the operator L in (|2.7p to a new 
operator of the same form 

L' = Lq — u', u'(x,x';q) = u'(x)S(x — x'), (3-1) 

given by means of an operator pair £, r\ according to the formulas 

L'C = (L, nL' = Lr,, (3.2) 

"twisting" L to L'. We consider the potential u(x) in L to be a real, smooth 
and rapidly decaying function of x and we look for self-conjugate £ and rj such 
that the new potential u'(x) is also real and smooth. In addition, we require r/ 
to be the left inverse of £, i.e., to obey the condition 

VC = I, (3.3) 

so that 

L = r,L'C (3.4) 

In order to get a new potential u'(x) not decaying along some directions of 
the plane, the two operator L and U must be related by a transformation more 
general than a similarity one. Therefore, we search for twisting operators having 
a product £77, which is not equal to /, but which is given by 

CV = I-P, (3-5) 

where P is an orthogonal self-conjugate projector, since as a consequence of (|3.3p 
P 2 = P and, thanks to the self-conjugate property of ( and 77, P* = P. The 
twisting operators £ and r\ generate a new potential v! via (13. 2)1 . but also the 
new dressing operators v 1 and ui'. In fact, taking into account (|2. 13)) . we get 
by d3HD L'Qv = C,Lv = (isL Q , and by (|2~T5)) ujrjL' = uoLri = L q lot]. Therefore, 
the operators v' and u>' defined by 

v' = (v, u' — cur), (3.6) 

obey equations 

LV = v'L , lu'L' = L lo', (3.7) 

analogous to (12.13)) satisfied by the dressing operators v and u. 

Note that because of ()3.3p the scalar product of these dressing operators is 
equal to / like for the original one in 1)2.14^ : 

wV = 1. (3.8) 

However, in contrast to v and lo, these operators do not satisfy a completeness 
relation, since, thanks to (|2 . 14[) and unlike it, one has by ()3.5p 

v'u' + P = I. (3.9) 



G 



In order to obtain a Darboux transformation we must specify the analyticity 
properties of the kernels q) and r](p; q) of the twisting operators with respect 
to the variables q. Since, thanks to (|2.14[) and (|3.6p , we have £ = v'lo and 
r] = uu)', the singularities of these kernels are related to the properties of v'(p; q) 
and u>'(p; q). We require that 

1. the kernels v'(p; q) and ui'{p; q) are independent of q 2 and have asymptotic 
behavior limq^oo u'(p; qj = 6(p) and limq^oo u'(p; q : ) = 5(p), 

2. the kernels f'(p;qi) and a/^q-J have (correspondingly, right and left) 
simple poles with respect to the variable q 1; i.e., there exist nontrivial 
limits 

v' h {p)= lim (qi -«6;)^'(p;qi), (3.10) 

q^ >ibi 

w'(p)= lim (qi+Pi -iajV(p;qi), (3.11) 

where ax, . . . , ajv a , &i, • ••, frjVj, are N a +Nb parameters, which we choose to 
be all different, and real in order to guarantee reality of the transformed 
potential, 

and, moreover, that 

3. the kernels ((p; q) and r)(p; q) besides the discontinuities at q 1Q = bi and 
q 1Q = Oj , which follow from requirement 2, have no additional departures 
from analyticity with respect to q x . 

Then, the kernels of the operators C, r\ are given by the following equations 

C(p; qi) 
»7(p; qi) 

where notations (|3 . 10[) and (|3.1ip were used. Let now %'(x,k) and £'(p;k) de- 
note the functions defined in terms of the dressing operators v' and u/ in analogy 
with (|2~T5|) and let x^fc) and £' aj (x) be their residua (cf. (|3~T0| and ([3~TTjl ). 
Then, it is easy to show that condition (|3.3|) is equivalent to the set of equations 

JV„ 

x!bi( x ) = -i^2x(z,iaj)™>ji(x), l = l,...,N b , (3.14) 

3=1 

£,0*0 = (x^foi&i), j=l,...,N a . (3.15) 



= s(p) dp 



1=1 



S(P) + Y, [ d P 



j v' hl {p-p')uj{p';ibi -p[) 

qi +P[ ~ ibi 
I v( p-p';ia j )uj' a .(p') 

qi +p'i - ia 3 



(3.12) 
(3.13) 



7 



where m(x) is the N a x A^-matrix with elements 

(a.j— bj)oo 

m j i{x)= I d yi e^- b ^-y^C.(y) X ' bl (y) , (3.16) 



which is well defined for bounded x' bl and £' a . and obeys 

d Xl mji(x) = (oj - bi)mji(x) - £' aj (x)x' bl {x). (3.17) 
In addition, from (|3.1[) and (|3.2[) . we get for the transformed potential u'(x) 

N a N b 

u'(x) = u(x) - 2d Xl ^2 X! £( x ' ib i)x( x > iaj)mji{x). (3.18) 

3=1 i=l 

Then, due to (|3.12[) and (|3 . 1 3[> . in order to define £ and 77 we must specify 
v' hi and uj,., i.e., X bl { x ) anc ^ £a( x )- Thanks to (|3.14[) and (|3. 15[> . this means 
that we have to define the matrix rriji(x). Inserting equations (|3 . 14[) and (|3.15[) 
into (|3.17p and using (|2.16[) , we write the following equation for the matrix 

fhji(x) = eW^-^'^mjtix) 

N a N b 

d Xl fhji(x) = - 2J '^mji l {x)^{x,ibv)^{x,ia r )fhj,i{x). (3.19) 

3'=U'=1 

Omitting details, we present, here, directly its solution in the two equivalent 
forms 

m(x) = {E Na + c.F(a;)) _:L c = c{E Nb + jF(x)c) _1 , (3.20) 

where c is an arbitrary real constant N a x N b matrix, Ejy a and E^ b are the 
unity N a x N a and N b x N b matrices, correspondingly, and J-(x) is a N b x N a 
matrix with elements !Fij{x) = !F(x,ibi,iaj), where 

xi 

F(x,k,k') = f dsi*(ar' J k)*(a',k') ; . (3.21) 

(k 3 -k' a )oo 

Now the potential (|3.18p also can be written in the two forms 

u'(x) = u(x) - 2dl x lndet(J5jv b + Pc) = u(x) - 2d 2 Xl \ndet{E Na + cT). (3.22) 

Inserting the matrix m found above into (|3.14j) . (|3 . 1 5[) and, then, the obtained 
equations into (|3.12[) . (|3.13p . we derive the explicit formulae for the dressing 
operators ( and rj. Finally, from (|3 .6[> using equation (|2 . 16[) and its analog for 
the transformed (primed) Jost and dual Jost solutions we get for them 

k) = k) - $(x, ia)(E Na + c J : (x))~ 1 c F(x, ib, k) = 

= ®(x, k) - $(x, ia)c{E Nb + Tix)^ 1 T{x, ib, k), (3.23) 
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f'(x,k) = $(x,k) - F(z,k,ia)(E Na + cT(x)y 1 c^(x,ib) = 

= V(x, k) - T(x, k, ia)c(E Nb + T{x) C y 1 ^{x, ib), (3.24) 

where j = 1, . . . , N a , 1 = 1,..., Nf,, and 

$(x,io) = diag{$(x,ia^)}, = diag{*(x,ify)}, 

.F(x, k, w) = diag{.F(x, k, «%■)}, .F(x, ib, k) = diag{.F(x, ib;, k)}, 

and where any of the two forms can be used. It is easy to see that both $'(x, k) 
and \&'(x, k) have poles at k = ibi and k = ia,-, correspondingly, and thanks 
to (|3.23p and (|3.24[) we get that the residua of these functions are given in terms 
of their values in the dual points by means of the relations 

^(x) = -ijVfoio,-)^, KM) =i£q,-,V(i6«), (3.25) 
i=i i=i 

as expected. One can show that the potential u'{x) in (|3 . 22[) and the Jost 
solutions from (|3 . 23[1 and (|3.24| coincide with those obtained in [11], including 
the case N a ^Nb considered here, which is recovered by choosing some zero rows 
or columns in the constant matrix C introduced in [11]. We postpone to Sec. [5] 
the discussion on the regularity conditions for the potential given by (|3.22[) . 



4 Resolvent 

Once the transformed operator L 1 has been obtained, the operators £ and rj can 
be used to investigate its spectral properties and existence and uniqueness of 
the corresponding resolvent M'. Multiplying the first and the second equation 
in (|3.2|) , respectively, from the right by 7/ and from the left by £, and recalling 
the definition (|3 - 5[) of P, we get the intertwining relation 

L' = (L V + L A , (4.1) 

where we introduced 

L A = L'P = PL'. (4.2) 
Thus we look for a resolvent in the form 

M' = (Mr) + M A , (4.3) 

where (Mr/ is determined by the above construction. Indeed, due to (|2.12l) 
and (|2.14[) we deduce that 

(Mr) = v'Mqw 1 , (4.4) 

and, then, we get for its kernel in the x-space the following bilinear expression 
in terms of the transformed Jost solutions 

{CMr)){x, x'-q) = - sgn(x 2 - 4)^—^ — / djuO ((g 2 +p\ - g?)(x 2 - 4)) x 
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x &{x; pi + iqi)f'(x'; p x + iq x ). (4.5) 

The second term in (|4.3I) . M\, is to be determined. By using (|3.2[) . (|2.9p . 
and (|3.5p , we deduce that M' in (|4.3|) is the right or left inverse of L' iff Ma is 
a solution, respectively, of the first or the second of the operator equations 

L'M A = P, M A L' = P. (4.6) 

Below we consider the solvability of these equations in the case of pure soliton 
potentials. In order to complete the discussion of the generic case, we mention 
that an explicit form of the operator P defined in (|3.5p can be derived by 
inserting there £ and 77 given in (|3 . 1 2[) and (|3 . 13[) . Thanks again to (|2.16|) . we 
get for the hat-kernel of this operator the expression 

M 

P(x,x';q)=i6(x 2 -x' 2 )y^e(qi-a n ) res $'(x,k)*'(x',k). (4.7) 
n=l 

Here we introduced the set of parameters 

{ai,...,a/s} = {ai,...,a Na ,bi...,b Nb }, M = N a + N b , (4.8) 

in order to make explicit the symmetry properties with respect to the parameters 
cij, bi used above. This kernel is different from zero only in the interval where q\ 
is between the lowest and highest values of the a m 's. Indeed, while for qi below 
this interval this is obvious by (|4.7[) . for the values of qi above the interval this 
follows from the equality 

N 

V res k)*'(a;', k) = 0, (4.9) 

n=l 

that in its turn is a consequence of (|3.25p . Now, according to the discussion 
above, we expect that Ma has a structure similar to P(x, x';q). It is clear that 
the kernel 

N 

M A (x,x';q) = Tid(±(x2-x' 2 ))Y9{ qi -a n ) res &(x, k)*'{x', k), (4.10) 

— ' k— ia n 

n=l 

obeys equations C' x M A (x, x'; q) = P(x,x';q), C' x , M A (x, x'; q) = P(x,x';q) for 
any sign. But it is easy to see that, in general, the kernel Ma(i,i';<;) con- 
structed by means of (|2.2[) is growing at space infinity and cannot be the kernel 
of an extended operator as defined in Sec. [2] Below we investigate this problem 
in detail by means of a specific example. 



5 Pure soliton potential and Jost solutions 

In the case u(x) = one gets the general iV-soliton solution, where N — 
maxima, Nb}. The transformed potential u'{x) and the corresponding Jost so- 
lutions can be easily obtained from the general expression derived in Sec. ([5]). 
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Thus, the transformed potential is given by any of the formulae in (|3.22[) where 
now the matrix F{x) has elements 



z (a j -b,)(x 1 + (a j +bi)x 2 ) 

an - bi 



(5.1) 



By expanding the determinants on the r.h.s. of (|3.22[) and by using the Binct- 
Cauchy formula one gets that, in order to have a regular potential, it is suffi- 
cient that the real matrix c satisfies the following characterization requirements 
(equivalent to those in [20]) 



A 



■ 7 

• ; Jn 



hi hi 



i Jn 



> 0, 



Ay = 



1 



k 



(5.2) 



for any 1 < n < minjA^, Nj,} and all minors, i.e., any choice of 1 < l\ < I2 < 
■ ■ ■ < l n < Nb and 1 < jt < j% < . . , j n < N a . Here, for the minors we used the 
standard notation 



3l,j2, ■ 

It, I2, ■ ■ 



7 l>n 



= dct 



a jlh 


a jih 


■ a h 




a,j 2 i 2 




a jnh 







(5.3) 



Formulae (|3.23p and (|3.24p , after rather cumbersome calculations that we skip 
here, can be more simply expressed as ratios of determinants. As in (|3.22[) one 
can alternatively, but equivalently, use determinants of (Nj, x Nf,)- or (N a x N a )- 
matrices. We use below the last case and we use the symmetric notation (|4.8[) . 
In these terms for the functions \' and £' related the Jost solutions by (|2. 16[) 
we get 




n , 6 , + ikr .W, 



(=1 

N h 



= U(b l+ ik) 



1=1 



t(x') ' 



where the t functions are determinants defined by 

r x (x,k) = det(%e A(x) (a + ik)T>), 
r e (x',k) = det^e^'^a + ik)- 1 ©) 
t(x) ^det(%e A ^ V), 



(5.4) 
(5.5) 



(5.6) 
(5.7) 
(5.8) 



i.e., determinants of (Nf, x A^-matrices obtained as products of square and 
rectangular matrices. Precisely, 



Bi, 



a 



N b -l 



a + i k = diag{a„ + i k} 



e* 4 ^ = diagie^^}, A n (x) = a n x x + 



a 2 n x 2 



(5.9) 
(5.10) 
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where i = 1,2, . . . , AT b , j = 1, 2, . . . , iV , rc. = 1, 2, . . . , M. As far as the potential 
is concerned, it can be expressed as u'(x) = — 2d* logr(x) and, taking into 
account (15. 8|) . one recovers the expression obtained in [19,20,22] using the tau 
functions. 



6 iV-soliton solutions in the case = 1 

Now, we restrict ourselves to the special case of the 7V-soliton potential with 
Nt, = 1 and N a = N arbitrary. In addition, for simplicity, without loss of 
generality, we choose a\ < a 2 < ■ ■ ■ < . In this situation the potential 
u'(x) has on the cc-plane N a ingoing rays and one outgoing ray, as schematically 
shown in Fig. 1. In this case r-functions (|5.6[) — (|5.8|) take the simple form 



■H * f P A(x) 

I] ;.,.-- 4(a:) (a m + ik), Tt(x,-k) = J2 h ^~-Z' 
z — ' * — ' n — + %k 



t x (x, k) = f m e^ x, (a m + ik), r ? (x, k) 

rn — 1 m— 1 



r(x) = f^ A{x) , (6-1) 



m— 1 



where the / m 's are obtained from the elements d m i's of matrix d in (|5. 1 1[) 
by a permutation which take into account the chosen ordering of the a m 's. 
Condition (|5.2[) here means that all f m > 0. 



Figure 1 : Rays and Polygon for N = 4 

We want to show that the extended operator L'{q) corresponding to this 
potential, as announced in the introduction, has a left self-conjugate annihilator 
K(q) for q belonging to a domain of the g-plane. Precisely, let us consider the 
domain inside the polygon inscribed in the parabola q 2 = q\ of the g-plane and 
with vertices at points q\ — a m , m — 1, . . . , N, whose characteristic function is 
given by 

rn — 1 
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x [9(0.2 - (a m + u n )qi + ot m a n ) - 8(q 2 - [a x + a^r)qx + a\aj^). (6.2) 

Notice that this polygon can be considered as dual to the ray structure (see 
Fig. 1 for the case M = 4). 

The main remark needed in order to get the annihilator K{q) is that the 
function 

n{q)e qx 

^(x\q) = — -r-— (6.3) 
t[x) 

is bounded in the x-plane. More precisely, as one can prove after a detailed 
study, it is exponentially decaying for x going to infinity in any direction on the 
plane, for q inside the polygon, and, for q on the borders of the polygon, it has 
directions of nondecaying (but bounded) behavior. 

Then, since thanks to (|2.16[) . (|5.4[) and (|6.ip 1/t is proportional to the value 
of the dual Jost solution ^'(a;,k) at k = ib\, we have 

£' x , d $(x';q)=0, (6.4) 

and, consequently, 

KL' = (6.5) 

for any operator K with kernel 

K(x, x'; q) = ip(x; q)ip(x';q) (6.6) 

where ip(x; q) is any arbitrary self-conjugate function, bounded in x and identi- 
cally zero outside the polygon introduced above. 

One can verify directly that KP = K . If, in addition, we choose ip(x; q) of 
the form ip(x; q) = (Pj)(x; q) with j(x; q) a self-conjugate function, bounded in 
x and identically zero outside the polygon, and such that 

dx~f(x;q)ip(x,q) = n(q), (6.7) 

one also get 

Ptp = ip, dx(p(x;q)i/](x,q) = n(q). (6.8) 



Then, it is easy to check that K is a self-conjugate projector commuting with 
the projector P, i.e., we have 

K* = K, K 2 = K, PK = KP = K. (6.9) 

The existence of this annihilator proves that the operator L' cannot have right 
inverse for q belonging to the polygon defined by the characteristic function (|6.2[) . 
In order to prove that, on the contrary, the resolvent M'(q) exists for q outside 
the polygon we go back to the hat-kernel Ma(i, x'; q) defined in (|4.10p and 
we study the boundedness properties of Ma(i, x'; q) = e~ q ( x ~ x ">Mi\(x,x']q). 
Since, thanks to (|4.9p . this function is identically equal to zero outside the strip 
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ai < <7i < OQV+i on the g-plane, we need to investigate its properties only in 
this strip. Using (|2. 16|) . (|5.4[) , (|5.5p . and (|6.1|) we rewrite it explicitly as 

±e-^-^)g(±(x 2 -4)) 

M A (a:,a; ;<j) = \ 

t(x)t(x') 

N 

x X! /™/^(<Zi-a™)(a™-«r l )e- 4 '" (x)+ - 4 " (a;) . (6.10) 

m,n=l 

Then, one can prove that, choosing the upper (bottom) sign for q above (below) 
the polygon in the strip ot\ < q\ < o/v+ij the kernel Ma (a;, x'; q) is a bounded 
function of x and q and it defines the kernel of an extended operator M&(q) 
according to the definition in Sec. El We conclude that, for q outside the polygon, 
Ma(<j) satisfy equations (|4.6|) and, consequently, the resolvent M'(q) exists and 
is given by (|4.3p . 

Notice that the total Green's function G(x,x',k) of the operator C! can 
be defined (see [11, 12]) as the value of the kernel M'(x,x';q) at q\ = kg, 
qi = — kpj. for a complex spectral parameter k = kj? +i kg . These values of 
q lie outside the parabola g 2 = <Zi, thus outside the polygon and touch it at the 
vertices only. Then the Green's function exists for any k but it is singular at 
the points k = ia m corresponding to the vertices of the polygon. These are the 
only singularities of the Green's function, since the discontinuities of the first 
term in ()4.3j) at qi = a m (see (|4.5p ) are compensated (outside the polygon) by 
the discontinuities of the second term, as follows from (|4.10j) . 

Finally, it is worth noting that the case N a = 1 and 7V& arbitrary can be 
handled in a analogous way, having the operator L'(q) a right instead of a left 
annihilator for q inside the polygon, and that in the case N a = Nb = 1 the 
polygon reduces to the segment of the line q 2 = («i + 02)91 — ai«2 with end 
points qi = ol\ and q\ = a 2 . In this case one recovers the results obtained 
in [11]. 
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